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Abstract 



In various models of super symmetry (SUSY), the lightest superparticle in the mini- 
^ \ mal SUSY standard model sector, which we call MSSM-LSP, becomes unstable. Then, 

h-^^ I we may observe the decay of the MSSM-LSP in the detector at the LHC experiment. 

^ ■ We show that the discovery of such a decay process (and the determination of the 

lifetime of the MSSM-LSP) may be possible at the LHC even if the decay length of 



t^^ I the MSSM-LSP is much longer than the size of the detector; sizable number of the 

^Q ■ MSSM-LSPs decay inside the detector if the lifetime is shorter than lO"*-^"^^ sec. We 

O ' also discuss the implications of the study of the MSSM-LSP decay for several well- 

motivated SUSY models. 

><" 



Supersymmetry (SUSY) is a well-motivated target of the LHC experiment. Indeed, not 
only the discovery but also detailed studies of the superparticles are possible at the LHC 
experiment if superparticles are within the kinematical reach [Tl[2j. In many of the studies, it 
is assumed that i?-parity is conserved, and that the lightest superparticle (LSP) is the lightest 
neutralino. If so, all the produced superparticles cascade down to the lightest neutralino just 
after the production and SUSY events are characterized by large missing px- 

Even though the lightest superparticle in the minimal-SUSY-standard-model (MSSM) 
sector, which we call MSSM-LSP, is often assumed to be the lightest neutralino and is 
stable, it is not always the case. Various scenarios where the MSSM-LSP becomes unstable 
have been proposed. One important example is the gauge-mediated SUSY breaking scenario 
[3], where the SUSY breaking in the MSSM sector is mediated via the standard- model 
gauge interaction so that SUSY-induced flavor violation are strongly suppressed. In the 
gauge-mediated model, the gravitino becomes the LSP and the MSSM-LSP decays into the 
gravitino. 

Another example is -R-parity violation, with which the MSSM-LSP becomes unstable. 
Usually, i?-parity conservation is assumed to realize the LSP dark matter scenario. However, 
LSP can be dark matter even if the i?-parity is violated; if the i?-parity violation is weak 
enough, the lifetime of the LSP becomes longer than the present age of the universe. This 
possibility becomes important when the gravitino is the LSP |4l |5] , because it has several 
advantages. In such a case, the thermal leptogenesis [^, which requires relatively high 
reheating temperature [7], may be possible without conflicting with the constraints from 
big-bang nucleosynthesis [8] and the overproduction of the gravitino [9]. In addition, if 
the gravitino is dark matter with i?-parity violation, a fraction of gravitino dark matter 
decays until the present epoch. The decay becomes a source of the high energy cosmic rays. 
In particular, recently, it has been shown that the anomalous excesses of the 7-ray and 
positron fluxes observed by EGRET [TU] and HEAT [llj experiments, respectively, can be 
simultaneously explained in the gravitino dark matter scenario with i?-parity violation if the 
lifetime of the gravitino is about 10^^ sec p^[T3] . In such a scenario, the MSSM-LSP decays 
mainly via the i^-parity violating interaction with the lifetime of 10"*^^"^^ sec. Discovery of 
the MSSM-LSP with such a lifetime may give us a hint to understand the origin of high 
energy cosmic rays. 

The experimental search for the decay of the MSSM-LSP gives important test of the 
scenarios with unstable MSSM-LSP. Since the superparticles are expected to be copiously 
produced at the LHC, we may have a chance to find the signal of the decay of MSSM- 
LSP. With long-lived unstable MSSM-LSP, we may see the decay at the LHC experiment 
in the form of a displaced vertex from the interaction point, non-pointing particle, and/or 
a disappearance of high pt tracks. Discovery of the MSSM-LSP decay is very important 
to understand the property of the MSSM-LSP. However, if the decay length of the MSSM- 
LSP is much longer than the size of the detector, most of the MSSM-LSPs escape from the 
detector before they decay. In such a case, the typical signal of the SUSY events are the 
same as the case where the MSSM-LSP is stable, and the discovery of the decay becomes 
statistically non-trivial. 
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In this letter, we discuss the possibihty of discovering the decay of the MSSM-LSP at 
the LHC experiment, paying particular attention to the case where the decay length is much 
longer than the size of the detector. We will show that the discovery may be possible if the 
lifetime is shorter than r ^ lO"*-^"^-' sec. Then, we consider the implication of the result for 
several types of models with unstable MSSM-LSP. Furthermore, if the decay is observed, 
the lifetime of the MSSM-LSP can be constrained p3l[T5] : we also discuss the possibility of 
determining the lifetime of long-lived MSSM-LSP. 

Let us start with discussing the basic formulae. If the MSSM-LSP is unstable, a fraction 
of the MSSM-LSPs produced at the LHC experiment decay inside the detector. If the 
MSSM-LSP (with its lifetime r) has the velocity v, the decay probability before propagating 
the distance L is given by Pdoc(-^) = 1 ~ q-l/v-tt ^ where 7 = (1 — v'^/c^)^^^'^ (with c ~ 
3.0 X 10^ m/sec being the speed of light). Then, denoting the pseudo-rapidity of the MSSM- 
LSP a.s rj = — lntan(^/2) (with 6 being the angle from the beam axis), the number of the 
MSSM-LSPs which decay inside the detector is given by 

iVdcc = iVtot / dr^dvfiv, v){l - e-^(— H.)Mr)^ (1) 



where iVtot is the total number of MSSM-LSP, l^^^'^^r]) is the distance to the outer boundary 
of the detector from the interaction point, and f{ri,v) = N^ldN/drjdv is the distribution 
function of the MSSM-LSP; J dr]dvf{r],v) = 1. 

In our following discussion, it is convenient to define 



L(^^) = I d7]dvf{7],v)'- 



Then, when the size of the detector is much smaller than the decay length ct, we obtain 

iVdec = iVtot . (3) 

CT 

When CT ^ NtotL'^'^^\ we expect several decay events inside the detector. Since, typically, 
f ~ c and 7 ~ 0(1 — 10), L*-^^^ becomes comparable to the size of the detector, as we will 
see in the following. Thus, roughly speaking, the number of the decay inside the detector is 
determined by the size of the detector and the total cross section for the SUSY events. In 
our analysis, for simplicity, we approximate the shape of the detector as a cylinder with the 
radius Ij^^^' and the half-length (to the 2r-direction) /(™'*^) = l)^^^' / tan 9 edge'- 

{,(max) I • n 
e-)/|cOS^| : V>Vedge ' ^ ^ 

where r^edge = — lntan(6'edge/2). From the muon chamber layout of the ATLAS detector, in 
our Monte Carlo (MC) analysis, we take [TJ^H 

4""^ = 10 m, r^edge = 1.0. (5) 



#^The end-cap of the ATLAS detector covers only up to jr/j < 2.7 jT|. We have checked that most of the 
MSSM-LSPs are within this region, and hence we do not impose a cut on rj for simphcity. 



Now, we calculate how many MSSM-LSPs decay inside the detector. The details depend 
on the MSSM parameters, and on what the MSSM-LSP is. When the MSSM-LSP is unstable, 
charged (or even colored) MSSM-LSP is phenomenologically viable. In the following, we 
discuss two of the important cases; one is the case where the lightest neutralino Xi is the 
MSSM-LSP while in the other case, the lighter stau f is the MSSM-LSP. For example, in 
the gauge-mediated model, they are two of the important candidates of the MSSM-LSP. 

First, we consider the case where the lightest neutralino Xi is the MSSM-LSP, and is a 
long-lived unstable particle. Even though Xi is invisible, spectacular signal may arise if Xi 
decays inside the detector. For example, in the models mentioned above (i.e., the gauge- 
mediated model or the model with i?-parity violation), Xi decays into 7 or Z boson and an 
invisible particle. If we can find a production of high energy 7 or the decay products of Z 
from the point which is displaced from the interaction point, like non-pointing photon [15] . 
it can be identified as the signal of the decay of Xi- (The study of these signals may require 
that the decay occurs at an inner region of the detector, which reduces the fiducial volume 
compared to Eqs. (jlj) and ([5]). For more details, see later discussion.) Signal of those events 
should be searched at the off-line analysis. In order to record those events, trigger may be 
an issue. One possibility is to use the missing p^ trigger. Because two XiS are produced in 
the SUSY events, even if one of Xi decays inside the detector to be the signal event, the 
other Xi is expected to escape from the detector. (Notice that we consider the case that ct 
is much longer than the size of the detector.) Such a non-decaying Xi should be a source 
of large missing pt- With the MC analysis, we calculate the distribution of the pt of Xi- 
The distribution for rrig = 1 TeV and 2 TeV are shown in Fig. [H With the present choice of 
parameters, more than 80 — 90 % of XiS have pt larger than 100 GeV for rrig = 1 — 2 TeV. 
Thus, assuming that the escaping x^ is the dominant source of the missing pt, most of the 
signal events have large missing pt- Thus, if a relevant missing pt trigger is implemented, 
the event can be recorded. 

In order to estimate how many Xi^ decay inside the detector, we perform MC analysis. At 
the LHC experiment, MSSM-LSP is mostly from the productions of colored superparticles: 
PP ~^ 99 1 M: aiid qq' (with g and q being the gluino and squark, respectively). The cross 
sections for these processes as well as the decay chains of the superparticles depend on the 
MSSM parameters. Here, as a well-motivated example, we adopt the gauge-mediated model 
to calculate these quantities. (Notice that, even though we use the gauge- mediated model 
to fix the underlying parameters, it is just for simplicity.) In our study, the SUSY events are 
generated as follows: 

1. Mass spectrum of the superparticles and their decay rates are calculated in the frame- 
work of the gauge-mediated model. Here, the simplest gauge-mediated model is 
adopted where the model is parametrized by tan/3 (i.e., the ratio of the vacuum ex- 
pectation values (VEVs) of up- and down- type Higgs bosons), number of vector-like 
messenger multiplets (in 5 + 5 representation of grand-unified SU{5) group) N5, the 
messenger scale Mmess? and A which is the ratio of the F-component of the SUSY 
breaking field to its VEV. 
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Figure 1: Distribution of pT of the final-state lightest neutralino with 10^ samples for the Xi' 
MSSM-LSP case. We take m-g = 1 TeV (shaded) and 2 TeV (unshaded). 

2. SUSY events expected at the LHC experiment are generated with y/s = 14 TeV. 

3. In each event, decay chains of superparticles are followed and the resultant momentum 
distribution of the MSSM-LSP is obtained. Then, L^^^^ is calculated. 

4. From the event samples, we calculate the number of the MSSM-LSP decay inside the 
detector as a function of r. 



In our analysis, we use ISAJET package [16] for the first step while, for second and third 
steps, HERWIG package [17] is used. 

For the study of the x?-MSSM-LSP case, we take tan/3 = 20, N^ = 1, M^^,, = 10^ GeV, 
and A is fixed to determine the mass scale of MSSM particles (in particular, in the present 
analysis, the gluino mass). With the above choice of parameters, the MSSM-LSP is Bino-like 
lightest neutralino, and its mass is given by m^o = 175 GeV, 360 GeV, and 500 GeV for the 
gluino mass rrig = 1 TeV, 1.5 TeV, and 2 TeV, respectively]^^! 

With the MC analysis, we found that L^^^^ is not sensitive to the mass spectrum of 
superparticles, and is ~ 10 m; L^^^^ = 8.5 m, 9.9 m, and 11.4 m for rrig = 1 TeV, 1.5 TeV, 



"^^Witli the present choice of parameters, the Higgs mass becomes smaller than the present experimental 

950 GeV. Because we choose the gauge-mediated model just as an 



bound of 114.4 GeV \m when 



< 



example of the SUSY model to fix the mass spectrum of MSSM particles, and also because the Higgs mass is 
sensitive to the masses of stops, we do not take the Higgs-mass constraint so seriously and extend our study 
to the parameter region of rUg ^ 950 GeV. If one is interested in the case of the gauge-mediated model, only 
the results for nig ^ 950 GeV are relevant. 
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Figure 2: The lifetime of Xi which gives A'dcc = 10 (with C = 100 fb ^) as a function of the gluino 
mass. 



and 2 TeV, respectively. We can see that L^°^^ is slightly enhanced as rrig increases, which 
is due to the decrease of the averaged velocity of the MSSM-LSP for larger value of rrig. On 
the contrary, the total cross section for the SUSY events strongly depends on the masses of 
superparticles; the total cross section is given by 1240 fb, 80 fb, and 10 fb, for rrig = 1 TeV, 
1.5 GeV, and 2 TeV, respectively. 

In Fig. [21 we plot the value of r which gives A'^dec = 10 (with C = 100 fb~^) as a function 
of the gluino mass. (Since A^dcc oc r^^, the lifetime which gives a different value of A'dec can 
be calculated from the figure.) For nig = 1 TeV, 1.5 TeV, and 2 TeV, A^dec > 10 requires r 
to be smaller than 7 x 10~^ sec, 6 x 10~^ sec, and 1 x 10~^ sec, respectively. Thus, when the 
lifetime is shorter than lO"^^"^-* sec, the number of the decay of the MSSM-LSP inside the 
detector turns out to be larger than ~ 0(1). 

Next, we consider the decay of f-MSSM-LSP. We expect several possibilities to find the 
signal of the decay of f. In the models we introduced, f decays into tau lepton and an 
invisible particle. For example, in the gauge-mediated model, the invisible particle is the 
gravitino while, in the scenario with the i?-parity violation, it is a neutrino. In those cases, 
we may find a displaced decay vertex by observing decay product (s) of the tau leptonj^^ 
In addition, if f propagates 0(10 cm) or so, hits in inner trackers should exist. Then, if 
we do not observe corresponding hits in calorimeters or in muon detector, we may identify 
such a short track as a signal of the f decay inside the detector. Furthermore, using the fact 
that two staus are produced, we can simply look for events with only one f track (which 



^^Otherwise, if selectron (smuon) is the MSSM-LSP, it decays into electron (muon) and an invisible particle. 
Then, energetic charged tracks from displaced vertices are the target. 
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Figure 3: The lifetime of f which gives A'^^ec = 10 as a function of the gluino mass (solid), and that 



which corresponds to N^^^ ™'"' 



10 with Vr, 



0.9 (dotted) and 0.8 (dashed). The integrated 



luminosity is taken to be £ = 100 fb 



may be useful to count the number of staus which decay very inside the detector for the 
determination of r). For the case of f-MSSM-LSP, the trigger may not be an issue because 
at least one high pt charged track (i.e., f) exists in the SUSY event; then, we expect that 
we can use the muon trigger. 

Here and hereafter, in the study of the case where f is the MSSM-LSP, we adopt the 
underlying parameters of tan/3 = 20, A^s = 2, and M^^ss = 10^ GeV, with which f becomes 
the MSSM-LSP. With such a choice of parameters, the mass of f becomes larger than 
100 GeV when m-g > 650 GeV, and is 170 GeV, 270 GeV, and 370 GeV, for m-g = 1 TeV, 
1.5 TeV, and 2 TeV, respectively. We follow the same procedure as the case of Xi-MSSM- 
LSP. For rrig = 1 TeV, 1.5 TeV, and 2 TeV, L^^^^ and the cross section for the SUSY events 
are given by 7.1 m and 2050 fb, 8.1 m and 150 fb, and 9.4 m and 20 fb, respectively. As 
in the case of Xi-MSSM-LSP, we can see that L^^^^ is insensitive to the mass spectrum of 
the superparticles. In Fig. [HI we show the lifetime which gives iVdec = 10. Again, when 



T 



< 



0(10 



-(3-5) 



sec), we can expect the decay of the MSSM-LSP inside the detector. 



Comparing Figs. [2] and [3l it is understood that iVdec does not depend much on the mass 
spectrum of the superparticles as far as the masses of the colored superparticles are fixed. 
This fact indicates that A^dec is mostly determined by the total cross section for the SUSY 
events (as well as by the lifetime of the MSSM-LSP). (See Eq. ([3]).) 

So far, we have seen that a sizable number of MSSM-LSPs decay inside the detector when 
r^ lO"'-'^"^-* sec. However, this does not necessarily mean that the decay of the MSSM-LSP 
can be easily observed. In particular, in our calculation of L^'^^\ we have assumed that the 
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typical size of the fiducial region is ~ 10 m. (More accurately, see Eqs. (jlj) and ([5]).) This 
implies that, in the case of the ATLAS detector, the decay of the MSSM-LSP may have to 
be identified by using only the last layer of the muon chamber. Such an analysis requires 
very high efficiency of detecting signals of the decay in the muon detector. In some cases, 
more conservative procedure to confirm the decay may be necessary. For example, for the 
case of the Xi-MSSM-LSP, some possibilities are (i) to require the decay of Xi inside the 
inner tracker region so that we can detect non-pointing photon using the electromagnetic 
calorimeter [T5], or (ii) to look for a decay in the calorimeters to see the decay vertex. For 
the case of f-MSSM-LSP, as we have mentioned, we can require (i) no hits in the muon 
detector for high pt charged particles observed by the inner detectors, or (ii) the discovery 
of the decay vertex of f in the calorimeters or inner trackers. Then, the size of the fiducial 
volume to study the decay is reduced. At the ATLAS detector, the inner radius of the 
muon chamber, the outer radius of the hadron calorimeter, and the outer radius of the inner 
tracker region are about 5 m, 4.25 m, and 103 cm, respectively [I]. Importantly, L^^^'^ given 
in Eq. ([2]) is proportional to the (typical) length of the fiducial region. Thus, if the MSSM- 
LSP is required to decay in the calorimeters or in the inner tracker, L^'^^^ is expected to be 
reduced by the factor of 2 — 10, and hence the maximum value of the lifetime with which 
sizable decay events are expected is. Even with such a smaller value of L^'^^\ as Figs. [2] and 
E] indicate, a significant number of the MSSM-LSP may decay in the fiducial region in the 
models we are interested in, as we see in the following. The understanding of the efficiency 
to find the signal of the decay should require extensive studies of the detector effects, which 
is beyond the scope of this letter. Here, we simply assume that the signals of the decay can 
be somehow identified in the following discussion. 

Let us now consider implications of the search for the decay of the long-lived MSSM-LSP. 
First, we consider the gauge-mediated model. When the Bino-like neutralino is the MSSM- 
LSP, it decays into gravitino and a gauge boson (i.e., photon or Z-boson), and the lifetime 
of Xi is estimated as 

where 77^3/2 is the gravitino mass. In the case where f is the MSSM-LSP, it decays into 
gravitino and the tau-lepton, and the lifetime of f is given by 



r 



2 X 10-^ sec X ( "^^ r' ( "^"^ y . (7) 

V200 GeVy VIOO keV/ ^ ' 



Thus, a sizable amount the MSSM-LSP decay is possible in both cases if m3/2 ^ 0(100 keV — 
1 MeV); in particular, using the fact that too large m^ is not preferred from the naturalness 
point of view, the number of the decay events can be as large as O(IOO) in such a case. 

Another important case is with i?-parity violation. In particular, with the introduction 
of the following bi-linear i?-parity breaking interaction into the SUSY breaking terms: 

>Crpv = ByipyLHu + h.c, (8) 



with L and Hu being the slepton and the up-type Higgs boson, respectively, it was pointed 
out that the EGRET and HEAT anomahes can be simultaneously explained if the gravitino 
is the LSP [T2l [13]. (Here and hereafter, we neglect the generation index for sleptons for 
simplicity.) Such a scenario works for any kind of the MSSM-LSP as far as the lifetime of the 
gravitino is ~ 10^^ sec. With the i?-parity violating interaction given in Eq. ([8]), gravitino 
%l)fj_ dominantly decays as t/^^ -^ vZ or l^W^ . When the gravitino is heavier than the weak 
bosons, the lifetime of the gravitino is estimated as [13] 



where k = Bbj>y/^^1 (with ma being the sneutrino mass) is the VEV of the sneutrino in units 
of the VEV of the standard-model like Higgs boson. With the i?-parity violation given in 
Eq. ([8]), the MSSM-LSP may also decay dominantly via the i?-parity violating interaction. 
Then, when the Bino-like neutralino is the MSSM-LSP, it decays into a neutrino and a 
standard-model boson with the lifetime 

r ~ 1 X 10"*^ sec x (-^) { ^^^} ,, ) ■ (10) 



,10-9/ V200 GeV. 

When f is the MSSM-LSP, it decays into the tau lepton and a neutrino, and the lifetime is 
given by 



"2 ^ 771^ \ ^ f rria ^2 



where m^ is the mass of the Bino-like neutralino. Requiring r3/2 ~ 10^^ sec to simultaneously 
explain the EGRET and HEAT anomalies, k is required to be ~ IQ-^ for 777,3/2 ~ 200 GeV, 
resulting in r ~ 0(10"^ sec) (0(10"^ sec)) when the MSSM-LSP is Bino-like neutralino 
(stau). From the study of A^dec given above, we expect sizable amount of the decay of 
the MSSM-LSP inside the detector with such a lifetime. In particular, when masses of 
superparticles are at most ~ 1 TeV to solve the naturalness problem, the number of the 
MSSM-LSP decay inside the detector is expected to be ~ 100 — 1000. Thus, the search 
for the decay of the long-lived MSSM-LSP should give an important test of the scenario to 
explain the origins of anomalous 7-ray and positron fluxes. 

So far, we have discussed the discovery of the decay of long-lived MSSM-LSP. Importantly, 
once the decay is found, we can also constrain the lifetime of the MSSM-LSP. This fact can 
be easily understood from Eq. ([3]); once the total number of the produced MSSM-LSP and 
the number of the decay in the detector are both determined, the lifetime r is constrained 
using Eq. ([3]). 

A relatively good determination of the lifetime is expected in particular when the MSSM- 
LSP is f (or other charged superparticles) even if the decay length is much longer than the 



size of the detector [H] |^ The tracking information about the long-lived charged particle 



'^^If the decay length ct of the MSSM-LSP is comparable or smaller than the size of the detector, decrease 
of the decay point may be experimentally observed. Then, the lifetime is determined by using the distribution 
of the decay points when ct^ a few cm |;19 . 



will be available if it travels transverse length longer than 0(10 cm). Thus, once large 
amounts of f samples become available, we expect that the distribution function / can be 
directly determined from the experimental data. Since we consider the case where ct is much 
longer than the size of the detector, most of the staus do not decay inside the detector, and 
are observed as energetic charged particles. Then, if enough amounts of fs are identified 
with the measurement of their momenta, we can determine /(r/, f). 

In discussing the identification of long-lived f tracks, it should be noted that the stau may 
be confused with the muon in particular when f does not decay inside the detector. This 
fact makes the determination of iVtot (or, more accurately, the number of events in which 
both of fs do not decay inside the detector) non-trivial. One possibility to distinguish f from 
the muon is to use the time-of-flight information; for this purpose, the transition radiation 
tracker and the muon system are useful in the ATLAS detector. If the velocity of f is small 
enough, it takes sizable time to reach trackers, calorimeters, and muon chamber after the 
production. If the velocity information is combined with the momentum information, f with 
small enough velocity can be distinguished from energetic muon. In our study, we require 
that the velocity of at least one of the two staus should be smaller than fmax so that the 
event can be identified as a SUSY event; we assume that, with the velocity measurement, f 
with V < fmax can be distinguished from high p^ muon whose velocity is almost the speed 
of light. We adopt several values of fmax to see how the results depend on fmax- On the 
contrary, if one of the staus decays inside the detector, it provides a striking signal, as we 
have discussed. 

Based on the above argument, we assume that the SUSY events can be identified if there 
exists at least one f which escapes from the detector with v < Vmax, or if one of fs decays 
inside the detector. (Since we consider the case ct ^ L^'^^\ we safely neglect the case 
where both of fs decay inside the detector.) Then, we expect that we can experimentally 
count the number of events where one of the f s decays inside the detector and the other 
escapes from the detector with its velocity smaller than fmax- We denote the number of such 
events as N^l^^^' . Importantly, N^lc^^"" can be calculated as a function of r after the 
experimental determination of the velocity distribution of f with the above-mentioned type 
of SUSY events. Thus, with the measurement of N^^ '^maxj^ ^ determination of r is possible. 
Since we are interested in the case where N^^^ '^"''"' is much smaller than the total number 
of f s observed, the statistical error in the determination of the lifetime is dominantly from 

If the correlation of the velocities of two staus is weak, N^^ ^™'''' is given by 

CT 

where N^^^ ""n^^^) jg ^j^g number of f s which satisfy v < Wmax, which is also experimentally 
measurable. With our MC analysis, we have confirmed that Eq. ( 1T2|) holds with a good 
accuracy. Thus, even though more accurate relation between -/V'dec^'''"'"' ^^^ r can be ob- 
tained once the experimental data become available, we use Eq. (IT^ to estimate how well 
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Figure 4: Distribution of v^jc of the final-state stau with 10^ samples. We take mg = 1 TeV 
(shaded) and 2 TeV (unshaded). 



the lifetime can be constrained for simplicity. 

In order to see the effect of the velocity cut, we show the distribution of v^ of f. The 
parameter A is taken to be 66 TeV and 141 TeV so that the gluino mass is given by vtig = 
1 TeV and 2 TeV, respectively. The result is shown in Fig. HI From the figure, we can see 
that f acquires large velocity in average. Thus, if one imposes relatively severe cut on the 
velocity of f to reduce the muon background, significant amounts of the stau events will be 
also discarded. (See the following discussion.) 



Next, we calculate A^tot ^^ 



0.8c and 0.9c; the lifetime which gives N^ot^'"'^^^' 



10 is also shown in Fig. [31 (Notice that A^dec corresponds to N^ 



(t)'<-(;max) 



dec 



with V„ 



Assuming that all the decay events are identified. Fig. [3] shows contours on rrig vs. r plane on 
which the statistical uncertainty in the determination of r is about 30 % (i.e., 5t = T/^/To). 
In order to determine r at this level, the lifetime is required to be, for fmax = c and 0.9c, for 
example, shorter than 1 x 10~^ and 4 x 10"^ sec (8 x 10~^ and 4 x 10~^ sec, 1 x 10~^ and 
8 X 10~^ sec) for rrig = 1 TeV (1.5 TeV, 2 TeV), respectively. Notice that, with a smaller 
input value of r, better determination of the hfetime is expected. 

As one can see, the sensitivity becomes worse as we require smaller value of fmax- This is 
because most of the staus are produced with relatively large velocity. (See Fig. HI) Impor- 
tantly, the muon system of the ATLAS detector is expected to have a good time resolution of 
about 0.7 ns, and the velocity measurement is possible with the accuracy of 6v/v'^ ^ 0.03c~^ 
[1]. Thus, good separation of f from /i may be possible even with a relatively large value 
of Vmax- In addition, as we have mentioned, the SUSY event may be identified even with- 
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out imposing the velocity cut if one of the staus decays inside the detector. Thus, if the 
total number of staus (with any velocity) can be somehow determined, the velocity cut is 
unnecessary. For example, if the correlation of the velocities of two staus is experimentally 
checked to be negligible, the total number of stau events can be calculated from the number 
of events with vi < fmax < ^2 and that of fi, f2 < fmax, where vi and V2 are velocities of two 
staus (with vi < V2). In such a case, we can adopt the result with v^aax = c. 

Such a measurement of the lifetime provide a quantitative information about underlying 
parameters. For example, in the gauge-mediated model, the measurement of the lifetime 
is possible when the gravitino is lighter than 0(100 keV — 1 MeV), and consequently, the 
gravitino mass can be determined assuming that f decays into the gravitino. Then, we 



can obtain an important information about the SUSY breaking scalej^ In addition, in the 
model with i?-parity violation, the size of the i?-parity violating coupling constant can be 
determined. 

Finally, we comment on the case where the lightest neutralino is the MSSM-LSP. Even 
though the discovery of the decay of (long-lived) x?; which is already very spectacular, may 
be possible, the precise determination of A^tot? which is necessary for the determination of 
the lifetime of Xi? is very non-trivial. This is because Xi is neutral and we cannot directly 
count the number of XiS which do not decay inside the detector. If the total number of 
SUSY events can be somehow estimated, however, the discovery of the decay of Xi gives 
significant information about the property of the MSSM-LSP. Once the SUSY events are 
found, we expect that the order-of-magnitude estimation of the number of SUSY events can 
be performed, which gives a bound on A^tot- Then, the measurement of A^dcc gives some 
information about the lifetime. 
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